Kinetic exchange opinion model: solution in the single parameter map limit 
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We study a recently proposed kinetic exchange opinion model (Lallouache et. al., Phys. Rev E 82, 
056112 (2010)) in the limit of a single parameter map. Although it does not include the essentially 
complex behavior of the multiagent version, it provides us with the insight regarding the choice of 
order parameter for the system as well as some of its other dynamical properties. We also study 
the generalized two-parameter version of the model, and provide the exact phase diagram. The 
universal behavior along this phase boundary in terms of the suitably defined order parameter is 
seen. 



I. INTRODUCTION 

Dynamics of opinion and subsequent emergence of con- 
sensus in a society are being extensively studied recently 
P-Q. Due to the involvement of many individuals, this 
type of dynamics in a society can be treated as an ex- 
ample of a complex system, thus enabling the use of con- 
ventional tools of statistical mechanics to model it 
Of course, it is not possible to capture all the diversities 
of human interaction through any model of this kind. 
But often it is our interest to find out the global per- 
spectives of a social dynamics, like average opinion of 
all the individuals regarding an issue, where the intrica- 
cies of the interactions, in some sense, are averaged out. 
This is similar to the approach of kinetic theory, where 
the individual atoms, although following a deterministic 
dynamics, are treated as randomly moving objects and 
the macroscopic behaviors of the whole system are rather 
accurately predicted. 

Indeed, there have been several attempts to realise the 
human interactions in terms of kinetic exchange of opin- 
ions between individuals [H, EH EH • Of course, there is no 
conservations in terms of opinion. Otherwise, this is sim- 
ilar to momentum exchange between the molecules of an 
ideal gas. These models were often studied using a finite 
confidence level, i.e., agents having opinions close to one 
another interacts. However, in a recently proposed model 
14], an unrestricted interactions between all the agents 
were considered. The single parameter in the model de- 
scribed the 'conviction' which is a measure of an agents 
tendency to retain his opinion and also to convince others 
to take his opinion. It was found that beyond a certain 
value of this 'conviction parameter' the 'society', made 
up of N such agents, reaches a consensus, where major- 
ity shares similar opinion. As the opinion values could 
take any values between [-1:+1], a consensus means a 
spontaneous breaking of a discrete symmetry. 



There have been subsequent studies to generalize this 
model, where the 'conviction parameter' and 'abilit y t o 
influence' were taken as independent parameters |15| . 
In that two-parameter version, similar phase transitions 
were observed. However, the critical behaviors in terms 
of the usual order parameter, the average opinion, were 
found to be non-universal. 

There have been other extensions in terms of a phase 
transition induced by negative interactions fl6| , an exact 
solution in a discrete limit 17], a generalized map ver- 
sion [lH and also the percolative transitions in a square 
lattices [l9j . 

In the present study, we investigate the single param- 
eter map version of the model, also proposed in Ref. [13] . 
In this limit, the map can be conceived as a random walk. 
Using standard random walk statistics, several static and 
dynamical quantities have been calculated. We show that 
the fraction of extreme opinion behaves like the actual 
order parameter for the system, and the average opinion 
shows unusual behavior near critical point. The critical 
behavior of the order parameter and its relaxation behav- 
ior near and at the critical point have been investigated 
both analytically and numerically. 



II. MODEL AND ITS MAP VERSION 

The kinetic exchange model of opinion pictures the opin- 
ion exchange between two agents like a scattering process 
in an ideal gas. However, unlike ideal gas, there is no 
conservation of the total opinion. This is similar to the 
kinetic exchange models of wealth distribution, where the 
conservation was also present [2lJ. The exchange equa- 
tions of the model read 



0*0 + 1) = XOi{t) + XeOj(t), 



(1) 
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and a similar equation for Oj(t + 1), where Oi(t) is the 
opinion value of the i-th agent at time t, A is the 'convic- 
tion parameter' (considered to be equal for all agents for 
simplicity) and e is an annealed random number drawn 
from a uniform and continuous distribution between [0:1], 
which is the probability with which i and j interact (see 
[HI). Note that the choice of i and j are unrestricted, 
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making the effective interaction range to be infinite. The 
opinion values allowed are bounded between the limits 
— 1 < Oi(t) < +1. So, whenever the opinion values 
are predicted to be greater (less) than +1 (-1) follow- 
ing Eq. (jTJ) , it is kept at +1 (-1). This bound, along with 
Eq. ([I]), defines the dynamics of the model. 

This model shows a symmetry breaking transition at 
a critical value of A (A c ~ 2/3). The critical behaviors 

N 

were studied using the average opinion O a — | X) Oi {t — * 

i=l 

oo)\/ N 12211 . An alternative parameter was also defined 
in Ref . [1J| , which is the fraction of agents having extreme 
opinion values. This quantity also showed critical behav- 
iors at the same transition point. 

The model in its original form is rather difficult to 
tackle a naly tically (it can be solved in some special limits 
though [13] )■ However, as it is a fully connected model, a 
mean field approach would lead to the following evolution 
equation for the single parameter opinion value (cf. [141 ]) 

0(t + 1) = A(l + e)0(t). (2) 

This is, in fact, a stochastic map with the bound \ 0(t)\ < 
1. For all subsequent discussions, whenever an explicit 
time dependence of a quantity is not mentioned, it de- 
notes the steady state value of that quantity and a sub- 
script a denotes the average over the randomness (i.e., 
ensemble average). As we will see from the subsequent 
discussions, this map can be conceived as a random walk 
with a reflecting boundary. As in the case of the multia- 
gent version, the distribution of e does not play any role 
in the critical behavior. We have considered two distri- 
butions, one is continuous in the interval [0:1] and the 
other is and 1 with equal probability. Both of these 
give similar critical behavior. 

We also briefly discuss the two-parameter model, 
where the 'conviction' of an agent and the ability to con- 
vince others were taken as two independent parameters 
[l5| . In that context, the map would read 

0(t + l) = {\ + pie)0{t), (3) 

where fi is the parameter determining an agent's ability 
to influence others. As before, |0(f)| < 1. 

III. RESULTS 
A. Random walk picture 

One can study the stochastic map described in Eq. J2j| by 
describing it in terms of random walks. Writing X(t) = 
\og(0(t)), Eq. ([2]) can be written as 

X(t + l)=X(t)+ V , (4) 

where, rj(t) = log[A(l + e)]. As is clear from the above 
equation, it actually describes a random walk with a re- 
flecting boundary at X — to take the upper cut-off of 
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FIG. 1: The average return time T of 0(t) to 1 in the map 
described in Eq. ([2| is plotted with (A — A c ). It shows a 
divergence with exponent 1 as is predicted from Eq. ([7]) 
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FIG. 2: O a is plotted with A. The data points are results 
of numerical simulations, which fits rather well with the solid 
line predicted from Eq. ([9]), with k — 0.7 

0(t) into account. Depending upon the value of A, the 
walk can be biased to either ways and is unbiased just 
at the critical point. As one can average independently 
over these additive terms in (|4l), this gives an easy way 
to estimate the critical point |20j. An unbiased random 
walk would imply (rj) = i.e., 

l 

y"log[A c (l + e)]de = (5) 
o 

giving A c = e/4, where we have considered an uniform 
distribution of e in the limit [0:1]. This estimate matches 
very well with numerical results of this and earlier works 
[l4| . In order to guess the A dependence of O a in the or- 
dered region, we first estimate the "average return time" 
T (return time is the time between two successive reflec- 
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tions) as a function of bias of the walk. For this uniform 

distribution of e, the average position to which the walker 

goes following a reflection from the barrier is (A + l)/2. 

The average amount of contribution in each step is given 
1 

by J" log [A (1 + e)]de — log(A/A c ). This, in fact, is a mea- 

o 

sure of the bias of the walk, which vanishes linearly with 
(A — A c ) as A — > A c . So, in this map picture, one would 
expect that on average by multiplying this A/A c factor T 
times, 0(t) would reach 1 from (A + l)/2. Therefore, 



A+l 



T~ =1, 



giving 



T 



log A 



log A - log A c 



log A 

A - A c 



(6) 



(7) 



for A — > A c . Clearly, the average return time diverges 
near the critical point obeying a power law: T ~ (A — 
A c ) -1 . In Fig. Q] we have plotted this average return time 
as a function of A. The power-law divergence agrees very 
well with the prediction. 

The steady state average value of X(t) i.e., X a (and 
correspondingly O a ) is expected to be proportional to 
VT in steps of log A: 

X a ~ VTlogX = kVTlogX, (8) 

where k is a constant. This gives 

O a = exp(-fc| logA| 3 / 2 (A - A c )- 1/2 ). (9) 

The above functional form fits quite well (see Fig. ^ 
with the numerical simulation results near the critical 
point. It may be noted that the numerical results for the 
kinetic opinion exchange (1) also fits quite well with this 
expression (9). We note that O a increases from zero at 
the critical point and eventually reaches 1 at A = 1. But 
its behavior close to critical point cannot be fitted with a 
power-law growth usually observed for order-parameters. 
Such peculiarity in the critical behavior of O a compels 
us to exclude it as an order parameter though it sat- 
isfies some other good qualities of an order parameter. 
Instead, we consider the average 'condensation fraction' 
p a as the order parameter. In the multi-agent version, 
it was defined as the fraction of agents having extreme 
opinion values i.e., -1 or +1. In this case it is defined as 
the probability that 0(t) = 1. We denote this quantity 
by p(t). As is clear from the definition, one must have 



Pa 



1 

T' 



(10) 



where, T is the return time of the walker. As T ~ (A — 
A c ) _1 ; Pa ~ (A— Ac)* 9 with /? = 1. This behavior is clearly 
seen in the numerical simulations (see Fig. 

Also, the relaxation time shows a divergence as the 
critical point is approached. We argue that there is a 
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FIG. 3: The average condensation fraction (probability that 
O = 1) p a is plotted with (A — Ac). A linear fit in the log-log 
scale gives the growth exponent 1, as predicted from Eq. (|10l) . 
Inset shows the variation of p with A. 



single time scale for both 0(t) and pit). So we calcu- 
late the divergence of relaxation time for 0(t) and nu- 
merically show that the results agree very well with the 
relaxation time divergence for both 0(t) and pit). Con- 
sider the subcritical regime where the random walker is 
biased away from the reflector and would have a proba- 
bility distribution as 

p[X) = -j, cxp(-B(X - vt) 2 )/t), (11) 

where w~ 1/T ~ (A— A c ) is the net bias and constants A, 
B do not depend on t. One can therefore obtain the prob- 
ability distribution P of O using p{X)dX = P(0)dO, 



P(0) = ^iexp(-S(logO 



vtf/t). (12) 



Hence 



O a it) 



OP{0)dO, 



A r 1 

-= \ exp(-BQogO-vt) 2 /t)dO 
yt Jo 



exp(-Sw 2 i), 

yt 



(13) 



in the long time limit, giving a time scale of relaxation 
r ~ v~ 2 ~ (A — A c ) -2 . We have fitted the relaxation of 
O a (t), obtained numerically, with an exponential decay 
and found r. As can be observed from Fig. [4] it shows a 
clear divergence close to critical point with exponent 2. 

We have obtained the relaxation time of p a {t) and it 
also shows similar divergence. Note that at A = A c , v — 
and it follows from Eq. that O a {t) - £ -1 / 2 . This be- 
havior is also confirmed numerically (see Fig. [S]). The 
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FIG. 4: The average relaxation time for O a {i) is plotted with 
A. This shows a prominent divergence as the critical point 
is approached. In the inset, the relaxation time is is plotted 
in the log-log scale against (A — A c ). The exponent is 2 as is 
expected from Eq. (|13[) 
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FIG. 5: The time dependence of both O a (t) and p(t) are 
plotted at the critical point in the log-log plot. The linear 
fit shows a time variation of the form t~ s with S = 1, as is 
expected from Eq. (fl3| . 

average condensation fraction p a (t) too follows this scal- 
ing, giving (5=1/2 (as order parameter relaxes as t~ s at 
critical point). 

We have also investigated the effect of having an ex- 
ternal field linearly coupled with 0(t). In the multiagent 
scenario, this can have the interpretation of the influence 
of media. The map equation now reads, 

0(t + l)=X(l + e t )0(t)+hO(t), (14) 

where h is the field (constant in time). We have studied 
the response of O a and p a at A = A c due to application 
of small h. We find that (see Fig. O both grows linearly 
with h. One expects the order parameter to scale with 
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FIG. 6: The variation of O a and p a are plotted against the 
external field h at the critical point A = Ac- The linear fit in 
log- log scale shows 5' = 1. 

external field at the critical point as p a ~ h 1 / 5 . In this 
case 8' = 1. 



B. Random walk with discrete step size 

One can simplify the random walk mentioned above and 
make it a random walk with discrete step sizes. This can 
be done by considering the distribution of e to be a double 
delta function, i.e., e = 1 or with equal probability. 
This will make t](t) in Eq. © to be log (A) or log(2A) with 
equal probability. Below critical point, both steps are in 
negative direction (away from reflector) and consequently 
taking the walker to — oo. Exactly at critical point (A = 
A c ) the step sizes become equal and opposite i.e., log A c = 
— log 2A C giving A c = \j\f2. Above critical point, one of 
the steps is positive and the other is negative. However, 
the magnitudes of the steps are different. This unbiased 
walker (probability of taking positive and negative steps 
are equal) with different step sizes can approximately be 
mapped to a biased walker with equal step size in both 
directions. To do that consider the probability p(x,t) 
that the walker is at position x at time t. One can then 
write the master equation 

p(x, t + 1) = ]-p(x + a,t) + ^p(x + a + b,t), (15) 

where a = log(A) and b = log(2). Clearly, 

dp(x,t) _ f b\ dp(x,t) fa 2 ab b 2 \ d 2 p{x,t) 
dt ~V + 2j dx + \Y + Y + JJ dx 2 

(16) 

One can compare this equation with the usual biased 
random walker equation 

dp(x,t) ,dp(x,t) {a') 2 d 2 p(x,t) 
—d^ = (P-^ a -d^ + - 2 ft?-" 
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FIG. 7: The comparison of the biased walk with equal step 
size with the original walk is shown. Reasonable agreement 
is seen for a wide range of A values. 



FIG. 8: The phase diagram for the two-parameter map as 
predicted from Eq. (I22|) (upper line). The phase diagram for 
the multi-agent version (fj, c = 2(1 — A c )) is also plotted for 
comparison (lower line). 



where p and q are the probabilities of taking a positive 
and negative steps (p + q = 1) and a' is the (equal) step 
size in either direction. Comparing these two equations 
one gets 



P = 




V(log(A)-log(A c ))2 + (log(A c ))2. (18) 



Therefore, as A — > A c , the bias {p — q) ~ (A — A c )/a'. 
These are consistent with the earlier calculations where 
we have taken the bias to be proportional to (A — A c ). To 
check if this mapping indeed works, we have simulated 
a biased random walk with above mentioned parameters 
and found it to agree with the original walk close to the 
critical point (see Fig. [7|). Similar to the approach taken 
for the continuous step-size walk, one can find the re- 
turn time of the walker. This time the walker is exactly 
located at A when it is reflected from the barrier. The 
return time again diverges as (A — A c ) _1 . Also O a (t) will 
have a similar form upto some prefactors. Condensation 
fraction will increase linearly with (A — A c ) close to the 
critical point. All the other exponents regarding the re- 
laxation time, time dependence at the critical point and 
dependence with external field are same as before. This 
shows that the critical behavior is universal with respect 
to changes in the distribution of e. 



of that two-parameter model. It would read 
0(t + l) = (X + ne)0(t). 



(19) 



As before, one can take log of both sides and in similar 
notations 



X(t+1) =X(t) + log(X + fj,e). 



(20) 



This can also be seen as a biased random walk. If one 
makes e discrete, this is again a walk with unequal step 
sizes, which can again be mapped to a biased walk with 
equal step sizes. Therefore, in either case, some pre- 
factors will be changed, but the critical behavior will be 
same as before. The critical behavior is, therefore, uni- 
versal when studied in terms of the condensation frac- 
tion. For uniformly distributed e (in the range [0:1]), one 
can get the expression for the phase boundary from the 
equation 



J log(A + fie)de = 0, (21) 



which gives 



log(A c + fi c ) H log — 

Me V A c 



(22) 



Of course, this gives back the logA c = e/4 limit when 
A = fi. The phase boundary is plotted in Fig. [5J It also 
agrees with numerical simulations. 



C. Two parameter map 



IV. SUMMARY AND CONCLUSIONS 



As the multi-agent model was generalized in a two- 
parameter model [l5j , one can also study the map version 



In this paper we have studied the simplified map ver- 
sion of a recently proposed opinion dynamics model. The 
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map can be cast in a random walk picture with reflecting 
boundary. Then using the standard random walk statis- 
tics, some steady state as well as dynamical behaviors 
are calculated and these are compared with the corre- 
sponding numerical results. It is observed that the usual 
parameter of the system, i.e. the average value of the 
opinion (O a ), does not follow any power-law scaling (see 
Eq. ©). In fact, it is the condensation fraction, or in this 
case the probability (p a ) that the opinion values touches 
the limiting value, turns out to be the proper order pa- 
rameter, showing a power-law scaling near critical point. 

The dynamical behaviors of the two quantities (O a {t) 
and p a it)) are similar. We have found the power-law re- 
laxation of these quantities at the critical point. Also, 



the divergence of the relaxation time on both sides of 
criticality shows similar behavior. We have also studied 
the effect of an external field in these models. At critical 
point, both the quantities grow linearly with the applied 
field. The average fluctuation in both of these quantities 
show a maximum near the critical point. These theoret- 
ical behaviour fits well with the numerical results. 

In summary, we develop an approximate mean field 
theory for the dynamical phase transition one observed 
for the map ([2]) and find the average condensation frac- 
tion p a ~ (A — \ C )P with (i — 1 behave as the order 
parameter for the transition and it has typical relaxation 
time t ^ (A — A c )~ 2 with z — 2 and at critical point 
A = A c (= e/4) decays as t~ s with 5=1/2. 
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